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We study nuclear effects in charged current deep inelastic neutrino-iron scattering in the frame- 
work of a analysis of parton distribution functions. We extract a set of iron PDFs and show 
that under reasonable assumptions it is possible to constrain the valence, light sea and strange 
quark distributions. Our iron PDFs are used to compute a;sj -dependent and Q^-dependent nuclear 
correction factors for iron structure functions which are required in global analyses of free nucleon 
PDFs. We compare our results with nuclear correction factors from neutrino-nucleus scattering 
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correction factors differ in both shape and magnitude from the correction factors of the models and 
charged-lepton scattering. 
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I. INTRODUCTION 

The high statistics measurements of neutrino deeply 
inelastic scattering (DIS) on heavy nuclear targets has 
generated significant interest in the literature since these 
measurements provide valuable information for global fits 
of parton distribution fimctions (PDFs) The use of 
nuclear targets is unavoidable due to the weak nature of 
the neutrino interactions, and this complicates the ex- 
traction of free nucleon PDFs because model-dependent 
corrections must be applied to the data. 

Additionally, these same data are also useful for 
extracting the nuclear parton distribution functions 
(NPDFs); for such an analysis, no nuclear correction fac- 
tors are required. Due to the limited statistics available 
for individual nuclear targets with a given atomic number 
A the standard approach is to model the ^-dependence 
of the fit parameters, and then combine the data sets 
for many different target materials in the global analysis 
0, S, i, S S, I3- However, the high statistics NuTeV 
neutrino-iron cross section data (> 2000 points) offer 
the possibility to investigate the viability of a dedicated 
determination of iron PDFs 

With this motivation, we will perform a fit to the 
NuTeV neutrino-iron data and extract the correspond- 
ing iron PDFs. Since we are studying iron alone and will 
not (at present) combine the data with measurements 
on different target materials, we need not make any as- 
sumptions about the nuclear corrections; this side-steps 
a number of difficulties 0, [13, [HI . 



While this approach has the advantage that we do not 
need to model the A-dcpendence, it has the drawback 
that the data from just one experiment will not be suffi- 
cient to constrain all the parton distributions. Therefore, 
other assumptions must enter the analysis. The theoret- 
ical framework will roughly follow the CTEQ6 analysis 
of free proton PDFs this will be discussed in Sec.HIl 

In Sec. mil we present the results of our analysis, and 
compare with nuclear PDFs from the literature. In 
Sec. IIVI we extract the nuclear correction factors from 
our iron PDFs and compare with a SLAC/NMC pa- 
rameterization taken from the £^-Fe DIS process 
and also with the parameterization by Kulagin & Petti 
fist ■ Finally, we summarize our results and conclu- 
sions in Sec. |Vl 



II. THEORETICAL FRAMEWORK 

A. Basic formalism 

For our PDF analysis, we will use the general features 
of the QCD-improved parton model and the analyses 
as outlined in Ref. Here, we will focus on the is- 

sues specific to our study of NuTeV neutrino-iron data 
in terms of nuclear parton distribution functions. We 
adopt the framework of the recent CTEQ6 analysis of 
proton PDFs where the input distributions at the scale 
Qo = 1-3 GcV are parameterized as [l^ 



i = Uy,d^,g,u + d,s,s . 
i = d/u , 



(1) 



r 



where u„ and dy are the up- and down-quark valence dis- 
tributions, u, d, s, s are the up, down, strange and anti- 
strange sea distributions, and g is the gluon. Further- 
more, the fi ~ /f denote parton distributions of bound 
protons in the nucleus A, and the variable < x < A 
is defined as a; :~ Ax a where xa = /'^PA ■ q is the 
usual Bjorken variable formed out of the four-momenta 
of the nucleus and the exchanged boson. Equation (jT]) is 
designed for < x < 1 and we here neglect^ the distribu- 
tions at a; > 1. Note that the condition fi{x > 1,Q) — 
is preserved by the DGLAP evolution and has the effect 

v/A 

that the evolution equations and sum rules for the f[ 



are the same as in the free proton case.^ 

The PDFs for a nucleus {A, Z) are constructed as 

ff{x,Q) = ^f!/^{x,Q)+^-^^fr^^{x,Q) (2) 

where wc relate the distributions inside a bound neutron, 
/"^^(a;,Q), to the ones in a proton by assuming isospin 
symmetry. Similarly, the nuclear structure functions are 



^ While the nuclear PDFs can be finite for x > 1, the magnitude 
of the PDFs in this region is negligible for the purposes of the 
present study (c/., Refs. 



^ While the quark number and momentum sum rules for the nu- 
clear case are satisfied as in the proton, there is no requirement 
that the momentum fractions carried by the PDF flavors be the 
same. A recent analysis at low found the Cornwall-Norton 
moments to be the same in iron as in deuterium formed from a 
free proton and a free neutron to within 3% Il6ll . 
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given by 



A 



{A-Z) 
A 



f:'^m) (3) 



such that they can be computed in next-to-leading order 
as convolutions of the nuclear PDFs with the conven- 
tional Wilson coefhcients, i.e., generically 



k 



(4) 



In order to take into account heavy quark mass effects we 
calculate the relevant structure functions in the ACOT 
scheme [13, [11 in NLO QCD [H. Finally, the differential 
cross section for charged current (anti-) neutrino-nucleus 
scattering is given in terms of three structure functions: 



G^ME 



dx dy 



Mxy i^,^^ 
2E 



^2xF,-^^±yil~^)xF,'' 



(5) 



It is instructive to compare this with the parton model 
expressions for the structure function F2 in l^A scatter- 
ing, where denotes a charged lepton: 



= ^iu^+u^) + l{d^+d^) + 



(14) 



Using the Callan-Gross relations in Eqs. ([7]) and (|12p . 
and neglecting the proton mass, the differential cross sec- 
tion Eq. ([5]) can be simplified in the form 

da (X (l-y + yV2)F/^±y(l-|)a;F/^ (15) 



with the limiting cases: 



da 



Ip V A I 1 p A 



(for y 



(for y^O) 



(16) 



where the ('— ') sign refers to neutrino (anti-neutrino) 
scattering and where G is the Fermi constant, M the 
nucleon mass, and E the energy of the incoming lepton 
(in the laboratory frame) . 



B. Constraints on PDFs 

We briefly discuss which combinations of PDFs can 
be constrained by the neutrino-iron data. For simplic- 
ity, we restrict ourselves to leading order, neglect heavy 
quark mass effects (as well as the associated production 
thresholds), and assume a diagonal Cabibbo-Kobayashi- 
Maskawa (CKM) matrix.^ The neutrino-iron structure 
functions are given by (suppressing the dependence on x 
and 



The latter form of da shows that the (anti-)neutrino 
cross section data naturally encodes information on the 

(-). (-). 

four structure function combinations i^2'^ ±xF^'^ and 



in separate regions of the phase space. 



If we assume s"" 



and 



the structure 



functions constrain the valence distributions d^ 



d^ - 



d^ 



7;A 



X 

1 



-F. 



uA 



u"' and the flavor-symmetric sea 
via the relations: 



2K 



Furthermore, we have 



(17) 
(18) 



Fi 

Fr 



uA 



jmjA 
2 

TPI^A 

^3 



d^' + s^ 
2xF'('^ , 
2[d^ + 



(6) 
(7) 

(8) 



The structure functions for anti-neutrino scattering are 
obtained by exchanging the quark and anti-quark PDFs 
in the corresponding neutrino structure functions: 



^1.2 - 

T^VA _ 



9] , 
9] • 



Explicitly this gives 



pi>A 

F^^ 
Fi^ 



2xF[ 
2 [i, 



vA 



(9) 
(10) 

(11) 
(12) 

(13) 



All these effects are properly included in our calculations. 



x 



(19) 
(20) 



Since we constrain the strange distribution utilizing the 
dimuon data,^ the latter two structure functions are use- 
ful to separately extract the d^ and d'^ distributions. 

For an isoscalar nucleus we encounter further simplifi- 
cations. In this case, u"^ 
implies — d^ — 



distributions are {v 



A f^A „A 



and ~ d"^ ~: which 
t 
}• 



V'". Hence, the independent quark 



It is 



Note that these equations are known not to be exact as the 
DGLAP evolution equations at NNLO generate an asymm etry 
even if one starts with s = s or c = c at some scale |20| . 
However, these effects are tiny and far beyond the accuracy of 
our study. 

5 See Refs. [H [H [H, [H, [H for details. 
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instructive to introduce the parameter A := 1/2 — ZjA 
which describes the degree of non-isoscalarity. This al- 
lows us to write the PDFs in a way which makes devia- 
tions from isoscalarity manifest: 



Ml, 



-dl'^] 



(21) 
(22) 
(23) 
(24) 



in terms of an averaged nuclear valence distribution 
= (u?^^ + (S'J^)I2 and an averaged nuclear sea dis- 
tribution = (vPl^ -\- (F/^)/2. Recall, fj"^ represents 
the distribution for a bound proton in the nucleus A] 
hence, the nuclear effects are encoded in these terms. 
Notice that non-isoscalar targets (A ^ 0) therefore pro- 
vide information on the difference between the valence 
distributions {u^^^ — d^^^) and the light quark sea dis- 
tribution (vP^^ — dP/^) in the nucleon. Unfortunately, 
the data are often corrected for non-isoscalar effects and 
this information is lost. 



Scheme 


Cuts 


Data 


# points 


2 

X 


xVpts 


Name 


ACOT 


Q > 1.3 GeV 


+ 


2691 


3678 


1.37 


A 




no Wcut 


V 


1459 


2139 


1.47 


kv 






V 


1232 


1430 


1.16 


kv 


ACOT 


Q > 2 GeV 


V \ V 


2310 


3111 


1.35 


A2 




W > 3.5 GeV 


V 


1258 


1783 


1.42 


A2i^ 






V 


1052 


1199 


1.14 


A2i^ 


MS 


Q > 1.3 GeV 


V -\-v 


2691 


3732 


1.39 


M 




no Wcut 


V 


1459 


2205 


1.51 








V 


1232 


1419 


1.15 




MS 


Q > 2 GeV 


V -\- V 


2310 


3080 


1.33 


M2 




W > 3.5 GeV 


V 


1258 


1817 


1.44 


M2i/ 






u 


1052 


1201 


1.14 


M2i^ 



TABLE I: Fits to NuTeV cross section and dimuon data. 



set of bands for fits of comparable quality (Ax^ ^ 50 for 
2691 data points) which provide an approximate measure 
of the constraining power of the data. 



C. Methodology 

The basic formalism described in the previous sections 
is implemented in a global PDF fitting package, hut with 
the difference that no nuclear corrections are applied to 
the analyzed data; hence, the resulting PDFs are for a 
bound proton in an iron nucleus. The parameterization 
of Eq. ID) provides enough flexibility to describe current 
data sets entering a global analysis of free nucleon PDFs; 
given that the nuclear modifications of the x-shape ap- 
pearing in this analysis are modest, this parameterization 
will also accommodate the iron PDFs. 

Because the neutrino data alone do not have the power 
to constrain all of the PDF components, we will need to 
impose some minimal set of external constraints. For ex- 
ample, our results are rather insensitive to the details of 
the gluon distribution with respect to both the overall 
and also the effect on the quark distributions. The 
nuclear gluon distribution is very weakly constrained by 
present data, and a gluon PDF with small nuclear modi- 
fications has been found in the NLO analysis of Ref. 0] ■ 
We have therefore fixed the gluon input parameters to 
their free nucleon values. For the same reasons the gluon 
is not sensitive to this analysis, fixing the gluon will have 
minimal effect on our results. Furthermore, we have set 
the dju ratio to the free nucleon result assuming that the 
nuclear modifications to the down and up sea are similar 
such that they cancel in the ratio. This assumption is 
supported by Fig. 6 in Ref. 

Because we have limited the data set to a single heavy 
target (iron), the ^ surface has some parameter direc- 
tions which are relatively flat. To fully characterize the 
parameter space, we perform many "sample fits" starting 
from different initial conditions, and iterate these fits in- 
cluding/excluding additional parameters. The result is a 



III. ANALYSIS OF IRON DATA 
A. Iron Data Sets 

We determine iron PDFs using the recent NuTeV 
differential neutrino (1371/1170 data points) and anti- 
neutrino (1146/966 data points) DIS cross section data 
Q where the quoted numbers of data points refer to the 
two different combinations of kinematic cuts introduced 
below. In addition, we include NuTeV/CCFR dimuon 
data (174 points) [2ll] which are sensitive to the strange 
quark content of the nucleon. 

There arc other measurements of neutrino-iron DIS 
available in the literature from the CCFR [13, [H, [H, [lO] , 
CDHSfH and CDHSW [13 collaborations; see, e.g., 
Ref. [33| for a review. There is also a wealth of charged 
lepton-iron DIS data including SLAC [H and EMC [sl, 
l36j.^ For the present study we limit our analysis to the 
NuTeV experiment alone; we will compare and contrast 
different experiments in a later study. 

B. Fit results 

The results of our fits to the NuTeV iron cross section 
and dimuon data are summarized in Table HI The cross 
section data have been corrected for QED radiation ef- 
fects, and the non-isoscalarity of the iron target [131; cor- 
respondingly, we have used A = 56, Z = 28 in Eqs. ([2]) 



^ C/. the Durham HEP Databases for a complete listing: 
http: / /www-spires . dur . ac .uk/hepdata/ 
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and ([3]) 7 Note, for an iron target the isoscalar correction 
factors are small and do not exceed the few % level. We 
have performed fits to the combined data as well as to 
the neutrino- and anti-neutrino data sets separately. Fur- 
thermore, two different cuts in the kinematic plane have 
been examined: a) Q > 1.3 GeV, no cut on the hadronic 
invariant mass W and b) Q > 2 GeV and W > 3.5 GeV, 
cf., Table m The NLO QCD calculation was performed 
in both the MS and ACOT schemes. The ACOT scheme 
calculation takes into account the heavy quark mass ef- 
fects, whereas the MS scheme assumes massless partons. 
The dominant target mass effects have been incorporated 
[3i,[3^.8 

As noted above, we have found bands for each class of 
fits from which we have chosen central representatives. 
The values have been determined taking into account 
the full correlations of the data employing the effective 
function given in Eq. (B.5) of Ref. |12l |. The numbers 
for the 'X^ /pts are roughly on the order of 1.4 for both the 
ACOT and the MS schemes.^ Furthermore, the fits to the 
anti-neutrino data have considerably better values; 
however, we will see below that this is at least partly due 
to the larger uncertainties of these data. 



CTEQ6.1M PDF's, the Base-2 fit does not apply any 
deuteron corrections to the data; this is in contrast to the 
Base-1 PDFs. Also, the Base-2 fit does include the CCFR 
data that has been corrected to a free nucleon using 
charged-lepton correction factors; the Fermilab CCFR 
experiment is the predecessor of NuTeV with comparable 
statistics as those from NuTeV HO]. The CCFR results 
in the large- a; region (x > 0.4) are consistently lower than 
those from NuTeV, and various sources contributing to 
the difference have been identified [1, [3] . One third of 
the discrepancy has been attributed to a mis-calibration 
of the magnetic field map of the muon spectrometer, i.e., 
to the muon energy scale in the CCFR analysis. About 
another third comes from model differences (cross sec- 
tion model, muon and hadron energy smearing models). 
A comparison of NuTeV and CCFR data can be found 
in Ref. Q. 

By comparing the free-proton PDF 'Base-1' and 
'Base-2' sets with the iron PDF sets of Table U we can 
gauge the size of the nuclear effects. Furthermore, differ- 
ences between observables using the 'Base-1' respectively 
the 'Base-2' reference sets will indicate the uncertainty 
due to the choice of the free-proton PDF.^^ 



1. PDF Reference Sets 

For the purposes of this study, we use two different ref- 
erence sets of free-proton PDFs which we denote 'Base-1' 
and 'Base-2'. 

Since we focus on iron PDFs and the associated nu- 
clear corrections, we need a base set of PDFs which are 
essentially free of any nuclear effects; this is the purpose 
of the Base-1 reference set [l^. Therefore, to extract the 
Base-1 PDFs we omit the CCFR and NuTeV data from 
our fit so that our base PDFs do not contain any large 
residual nuclear corrections.^" The absence of such nu- 
clear effects will be important in Sec lIVI when we extract 
the nuclear corrections factors. 

The Basc-2 PDFs arc essentially the CTEQ6.1M PDFs 
with a modified strange PDF introduced to accommo- 
date the NuTeV dimuon data.^^ In the manner of the 



We have checked that omitting the isoscalar correction factors 
and using A = 56, Z = 26 gives almost identical results. 

* Target mass effects (TMC) are expected to be relevant at large 
Bjorken-x or small momentum transfers |39||. For issues of 
higher orders and higher twist cf. Refs. [40l. |41|. 13 l43l| 

® Fits to this same data neglecting the correlations between the 
errors and using the conventional function {cf. Eq. (B.l) in 
[12|| ). have smaller /pts ~ 1. While the uncorrelated errors are 
larger, the extracted parameters are similar. 

" We do retain the deuteron data as this has onlya small correction 
over the central x-range, {cf Sec. UVA)! g, llOl . The deuteron 
correction has been applied in the Base-1 fit. Also, for the DrcU- 
Yan Cu data (E605), the expected nuclear corrections in this 
kinematic range are small (a few percent) compared to the overall 
normalization uncertainty (15%) and systematic error (10%). 

^ These PDFs have been determined from a fit to the same data 



2. Comparison of the Fits with Data 

The quality of our fits of Table U can also be observed 
directly in Figures [1] - [3] where we compare the theoret- 
ical cross section {\ / E)(P a / dx dy with a selection of the 
data. To be specific, we show all the data taken with 
beam energies E = 65, 150, and 245 GeV which pass our 
kinematic cuts. The measurements are organized in bins 
of function of the inelasticity y and cover the a;- 

range 0.015 < x < 0.750. The momentum transfers can 
be computed using the relation = 2MExy. We nor- 
malize these plots using the 'A2' fit which implements 
the kinematic cuts Q > 2 GeV and W > 3.5 GeV (c/. 
Table m. We note that these are the cuts employed in 
the CTEQ6 analysis in order to reduce the sensitivity to 
target mass and higher twist effects. 

The fit provides a good description of the data which 
are distributed around unity for most of the bins. For 
reference, the results of fit 'A' (solid line) and Base-1 



set as in the CTEQ6 analysis with the addition of the the NuTeV 
dimuon data. The changes to the strange sea induce only minor 
changes to the other fit parameters; this has a minimal effect on 
the particular observables {da, F2) we examine in the present 
study. 

All results have been computed with both Basc-1 and Base-2 
PDFs. Since the Base-2 PDFs use CCFR and NuTeV data, the 
resulting PDFs will depend on the nuclear corrections which we 
are trying to determine. Therefore, we will predominantly display 
the Base-1 PDFs for comparison in the following Sections. 
Conversely, global analyses of nuclear PDFs tend to use looser 
kinematic cuts due to the lack of small-x data and the interest 
in the very large- a; region. 
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FIG. 1: Representative comparison of fit 'A2' to the NuTeV neutrino and anti-neutrino cross section data. Shown are the 
data points for various a;-bins versus the inelasticity y for an energy of iJ = 65 GeV in a data-over-theory representation. For 
comparison, we also show results for the Base-1 PDFs (dotted) and the 'A' fit (solid); the fit 'A2' imposes more stringent cuts 
on Q > 2 GeV and W > 3.5 GeV. 



PDFs (dotted line) are shown as well. For fit 'A2', the 
effect of the Q > 2 GeV cut is to remove data at low 
y in the small-x region, and the > 3.5 GeV cut ex- 
cludes \ow-y data at large x. The effects of these cuts 
on the fit are visible by comparing the difference of the 



sohd hne ('A') from unity ('A2'). For x > 0.045, we ob- 
serve minimal differences between the 'A' and 'A2' fits, 
and conclude the effect of the kinematic cuts (Q > 2 GeV 
and W > 3.5 GeV) arc nominal in this region. In the low- 
est X bin {x ~ 0.015), much of the data is eliminated by 
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FIG. 2: The same as in Fig. [T]for a neutrino energy of E = 150 GeV. 



the Q > 2 GeV cut such that fit 'A2' is only constrained 
by a few data points at large y for the higher neutrino 
energies, cf. Fig. [3l Since both, fit 'A' and fit 'A2', have 
large uncertainties in this x-region the comparison of in- 
dividual representatives is less significant — in particular 
at medium and low y where no data points lie. In con- 
clusion, we discern no relevant differences between the 
two classes of fits over the entire kinematic plane and 



will therefore mainly focus on fit 'A2' in the following 
sections. 



3. Comparison of the Fits with Reference PDFs 

The dotted curve in Figures [T] - [3] shows the cross sec- 
tions obtained with Base-1 free-proton PDFs, inserted 
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FIG. 3: The same as in Fig. [T]for a neutrino energy of E = 245 GeV. 



into Eq. ^ to obtain "free iron" PDFs, divided by the 
cross sections computed with fit 'A2' PDFs. The Base-2 
PDFs (not shown) yield similar results as we demonstrate 
in Sec. IIV CI We expect the base PDFs will provide a 
poorer description of the data since the nuclear modifi- 
cations are not taken into account; the deviations of these 
curves from unity indicate the size of the nuclear effects. 

We observe that the Base-1 results at small-x {x ~ 



[0.045 - 0.08]) are generally below unity (the 'A2' fit) in 
the y region of the data points implying an enhancement 
due to nuclear effects. As discussed above, the results in 
the lowest x bin {x = 0.015) are less clear as the uncer- 
tainties are larger since the kinematic cuts remove much 
of the data. Nevertheless, do not see a clear signal of 
shadowing in this region (c/., Fig. [3] at large y). 

For intermediate x ~ [0.125 — 0.175] the Basc-1 (dotted 
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line) results are very similar to fit 'A2'. For larger x ~ 
[0.225 — 0.65] we observe a suppression of the nuclear 
cross sections qualitatively similar to what is known from 
charged lepton DIS. Finally, in the region a; > 0.75 the 
nuclear cross section is again enhanced — an effect usually 
attributed to the Fermi motion of the nuclcons in the 
nucleus. 

In conclusion, we observe the following pattern for the 
nuclear cross section compared to the free nucleon cross 
section: i) enhancement for x > 0.75, ii) suppression for 
X - [0.225 - 0.65], iii) equality for x - 0.125, and iv) 
slight enhancement for x ~ [0.045 — 0.08]. This is to 
be contrasted with the expectation from charged lepton 
DIS with the well-known pattern: i) enhancement for x > 
0.75 (Fermi motion), ii) suppression for x ~ [0.3 — 0.8] 
(EMC effect), iii) enhancement for x ~ [0.06-0.3] (Anti- 
shadowing); and iv) suppression for x < 0.06 (Shad- 
owing). Thus, for X > 0.3 our results arc generally as 
expected. However, wc find that the usual behavior at 
medium and small x is modified. We will examine this 
further in the following sections. 



C. Iron PDFs 

Having established the quality of our fits, we now ex- 
amine the nuclear (iron) parton distributions ff{x,Q'^) 
according to Eq. Figure 2] shows the PDFs from 

fit 'A2' at our input scale Qq = rric = 1.3 GeV versus 
X. For an almost isoscalar nucleus like iron the u and d 
distributions are very similar, see Eqs. ((2T|) - p4)) . There- 
fore, we only show the Uy and u partons, together with 
the strange sea.^^ As explained above, the gluon dis- 
tribution is very similar to the familiar CTEQ6M gluon 
at the input scale such that wc don't show it here. In 
order to indicate the constraining power of the NuTcV 
data, the band of reasonable fits is depicted. The fits in 
this band were obtained (as outlined above) by varying 
the initial conditions and the number of free parameters 
to fully explore the solution space. All the fits shown 
in the band have / DOF within 0.02, which roughly 
corresponds to a range of Ax^ ~ 50 for the 2691 data 
points. 

As can be seen in FigurelH the distribution (Fig.[4|i) 
has a very narrow band across the entire x-range. The 
up- and strange-sea distributions (Fig. [4}3 and Fig. ^) 
arc less precisely determined. At values of x down to, 
say, X ~ 0.07 the bands are still reasonably well con- 
fined; however, they open up widely in the small- a; re- 
gion. Cases where the strange quark sea lies above the 
up-quark sea are unrealistic, but are present in some of 
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FIG. 4: Parton distributions for iron at our input scale 
= 1.69 GeV^ Shown are the bands (in yellow) from fit 
'A2' for the up quark valence distribution (upper figure), the 
up quark sea (middle), and the strange quark sea (lower fig- 
ure). The central PDF from fit 'A2' is shown by the solid line. 
The dashed lines depict parton distributions constructed ac- 
cording to Eq. ([2]) with j4 = 56 and Z = 2iS using the Base-1 
free-proton PDFs. The dotted lines are the leading order 
HKN04 nuclear parton distributions [3], the dotted-dashed 
lines are the next-to-leading order (NLO) HKN07 nuclear par- 
ton distributions Q], and the dot-dashed lines are the next- 
to- leading order distributions (DS) from Ref. 0. The vertical 
line marks the lower limit of the data in the x variable. 



While iron is roughly isoscalar, other nuclear PDFs can exhibit 
larger differences between the u and d distributions — the ex- 
treme case being the free-proton PDF. When comparing PDFs of 
Eq. JJll, we must keep in mind that it is ultimately the structure 
functions defined by Eq. l[4]l which arc the physical observables. 



the fits since this region {x < 0.02) is not constrained by 
data. We have included the curves for our free-proton 
Base-1 PDFs (dashed), as well as the HKN04 Q (dot- 
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ted), the NLO HKN07 Q (dotted-dashed) , and DS 
(dot-dashed) nuclear PDFs.^^ 

The comparison with the Base-1 PDFs is straightfor- 
ward since the same theoretical framework (input scale, 
functional form, NLO evolution) has been utilized for 
their determination. Therefore, the differences between 
the solid band and the dashed line exhibit the nuclear 
effects, keeping in mind that the free-proton PDFs them- 
selves have uncertainties. 

For the comparison with the HKN04 distributions, it 
should be noted that a SU(3)-flavor symmetric sea has 
been used; therefore, the HKN04 strange quark distri- 
bution is larger, and the light quark sea smaller, than 
their Base-1 PDF counterparts over a wide range in x. 
Furthermore, the HKN04 PDFs are evolved at leading 
order. 

In a recent analysis, the HKN group has published a 
new set of NPDFs (HKN07) including uncertainties [1]. 
They provide both LO and NLO sets of PDFs, and we 
display the NLO set. These PDFs also use a more general 
set of sea distributions such that u{x) ^ d{x) ^ s{x) in 
general. 

The DS PDFs are hnked to the GRV98 PDFs [H] with 
a rather small radiatively generated strange sea distri- 
bution. Consequently, the light quark sea is enhanced 
compared to the other sets. Additionally, the DS sets 
are evolved in a 3-fixed-flavor scheme in which no charm 
parton is included in the evolution. However, at the scale 
Q = nic of Fig. [4] this is of no importance. 



IV. NUCLEAR CORRECTION FACTORS 

In the previous section we analyzed charged current 
I'-Fe data with the goal of extracting the iron nuclear 
parton distribution functions. In this section, we now 
compare our iron PDFs with the free-proton PDFs (ap- 
propriately scaled) to infer the proper heavy target cor- 
rection which should be applied to relate these quantities. 

Within the parton model, a nuclear correction factor 
R[0] for an observable O can be defined as follows: 



R[0] 



^[NPDF] 
0[free] 



(25) 



where ©[NPDF] represents the observable computed 
with nuclear PDFs, and ©[free] is the same observable 
constructed out of the free nucleon PDFs according to 
Eq. (P5|) . Clearly, R can depend on the observable under 
consideration simply because different observables may 
be sensitive to different combinations of PDFs. 



This means that the nuclear correction factor R for 
and will, in general, be different. Additionally, the 
nuclear correction factor for will yield different results 
for the charged current v-Fe process {W^ exchange) as 
compared with the neutral current £^-Fe process (7 ex- 
change). Schematically, we can write the nuclear cor- 
rection for the DIS structure function F2 in a charged 
current (CC) v-A process as (cf. Eq. d?])):^^ 



R''cc{F2;x,Q'') 



(26) 



and contrast this with the neutral current (NC) £^-A 
process (cf. Eq. (fH|l ): 

R%''ciF2;x,Q^)^ 



(27) 



where the superscript " 0" denotes the "free nucleon" 
PDF which is constructed via the relation: 

ffix,Q) = jfrix,Q)+^-^-j^frix,Q) . (28) 

Clearly, the i?-factors depend on both the kinematic vari- 
ables and the factorization scale. Finally, we note that 
Eq. (|25[) is subject to uncertainties of both the numerator 
and the denominator. 

We will now evaluate the nuclear correction factors for 
our extracted PDFs, and compare these with selected 
results from the literature [l^, [13, [iBl-^^ Because we 
have extracted the iron PDFs from only iron data, we 
do not assume any particular form for the nuclear A- 
dependence; hence the extracted R[0] ratio is essentially 
model independent. 



A. Deuteron corrections for the F^f'^/F^ ratio 

The structure function ratio provides a com- 

mon (and useful) observable to use to gauge the nuclear 
effects of iron. To construct the numerator, we will use 
our iron PDFs as extracted in fits 'A' and 'A2.' For 
the denominator, we will use the Base-1 and Base-2 free 
proton PDF; however, converting from free proton struc- 
ture functions to deuteron structure functions is nontriv- 
ial and model-dependent. 

In Fig. [5] we display the NMC data for F^/F^ ^ 
and compare this to a variety of data parameterizations 



In a recent publication, Eskola et al. perform a global reanal- 
ysis of their ESK98 nuclear PDFs. While we do not present a 
comparison here, the results are compatible with those distribu- 
tions displayed in Fig. |4] a comparison can be found in Figs. 10 
and 11 of Ref. 0. 



^° The corresponding anti-neutrino process is obtained with a u <-» 
d interchange. 

Note that our comparison with the Kulagin-Petti model is based 
on the work in Ref. Il4l . 
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FIG. 5: NMC data for F^/F^ ^ a.t = 5.47 GeV^ in 
comparison with the theory prediction for F^^ /F2 computed 
using free-proton Base-2 PDFs. The dashed line shows the 
structure function ratio obtained with the Base-1 PDFs; in 
this case a nuclear correction factor for deuterium has been 
applied (c/., Refs. 

For comparison, we also show the 
parameterizations of Arneodo et al. [4^ and Tvaskis et al. 



0, m, 113, Si- The dashed line shows the structure 
function ratio computed with the Base-1 PDFs; in this 
case a nuclear correction factor for deuterium has been 
applied using the parameterization of Ref. [l^l . The solid 
line shows the structure function ratio computed with the 
Base-2 PDFs; in this case no nuclear correction factor for 
deuterium was applied. The dotted line (Arneodo) is the 
parameterization of Ref. [i^, and the dot-dashed line 
(Tvaskis) is the parameterization of Ref. [13, We 
see that the range of discrepancies in the deuterium cor- 
rections are typically on the order of a percent or two 
except at large x; while this correction cannot be ne- 
glected, it is small compared to the much larger iron 
nuclear corrections. To explore a range of possibilities 
(reflecting the underlying uncertainty) we have incorpo- 
rated deuteron corrections into the Base-1 PDF, but not 
the Base-2 PDF; hence the spread between these two ref- 
erence PDFs will, in part, reflect our ignorance of and 
other uncertainties of proton PDFs at large- a;. 



B. 



Fp 



IF? 



for neutral current (NC) charged 
lepton scattering 



We will also find it instructive to compare our re- 
sults with the F^^jF? as extracted in neutral current 
charged-lepton scattering, l^-Fe. In Fig. [6] we compare 
the experimental results for the structure function ra- 
tio F^^jF^ for the following experiments: BCDMS-85 
g^, BCDMS-87 SLAC-E049 [H, SLAC-E139 [|, 
SLAC-140 [H]. The curve (labeled SLAC/NMC param- 
eterization) is a fit to this data [l^. While there is a 
spread in the individual data points, the parameteriza- 
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FIG. 6: Parameterization for the neutral current charged lep- 
ton structure function F^^ jF^ ■ For comparison we show ex- 
perimental results from the BCDMS collaboration (BCDMS- 
85 [111, BCDMS-87 H) and from experiments at SLAG 
(SLAC-E049 m, SLAC-E139 [|, and SLAC-E140 [3). Nor- 
malization uncertainties of the data have not been included. 



tion describes the bulk of the data at the level of a few 
percent or better. It is important to note that this pa- 
rameterization is independent of atomic number A and 
the energy scale [S^]; this is in contrast to the results 
we will derive using the PDFs extracted from the nuclear 
data.^^ Additionally, we note that while this parameter- 
ization has been extracted using ratios of F2 structure 
functions, it is often applied to other observables such 
Fi^3^L or da. We can use this parameterization as a 
guide to judge the approximate correspondence between 
this neutral current (NC) charged lepton DIS data and 
our charged current (CC) neutrino DIS data. 



C. Correction Factors for d^a/dxdQ^ 

We begin by computing the nuclear correction factor 
R according to Eq. ((25| for the neutrino differential cross 
section in Eq. ^ as this represents the bulk of the NuTeV 
data included in our fit. More precisely, we show R- 
factors for the charged current cross sections dPa/dx dQ^ 
at fixed which can be obtained from Eq. ([S]) by a 
simple Jacobian transformation and we consider an iron 
target which has been corrected for the neutron excess, 
i.e., we use the PDFs in Eq. ^ (for the numerator) and 
Eq. dlH) (for the denominator) with A = 56 and 2" = 28. 
Our results are displayed in Fig. [7] for = 5 GeV^ 
and a neutrino energy Ei, = 150 GeV which implies, 
due to the relation = 2MEj^xy, a minimal x-value 



In particular, we will find for large x (> 0.5) and Q comparable 
to the proton mass the target mass corrections for /F^ are 
essential for reproducing the features of the data; hence the Q 
dependence plays a fundamental role. 
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FIG. 7: Nuclear correction factor R according to Eq. (|25[) for 
the differential cross section cPa/dxdQ^ in charged current 
vFe scattering at = 5 GeV^ and Ev = 150 GeV. Results 
are shown using the 'A2' fit for the charged current neutrino 
(solid lines) and anti-neutrino (dashed lines) scattering from 
iron. The upper (lower) pair of curves shows the result of 
our analysis with the Base-2 (Base-1) free-proton PDFs. The 
correction factors shown here are for an iron target which has 
been corrected for the neutron excess. 



of Xniin = 0.018. The solid (dashed) lines correspond to 
neutrino (anti-neutrino) scattering using the iron PDFs 
from the 'A2' fit. 

We have computed R using both the Base-1 and Base-2 
PDFs for the denominator of Eq. ((25|) ; recall that Base- 
1 includes a deuteron correction while Base-2 uses the 
CCFR data and does not include a deuteron correction. 
The difference between the Base-1 and Base-2 curves is 
approximately 2% at small x and grows to 5% at larger x, 
with Base-2 above the Base-1 results. As this behavior is 
typical, in the following plots (Figs. [Hand Figs.[9l) we will 
only show the Base-1 results. We also observe that the 
neutrino (anti-neutrino) results coincide in the region of 
large x where the valence PDFs are dominant, but differ 
by a few percent at small x due to the differing strange 
and charm distributions. 
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FIG. 8: Nuclear correction factor R according to Eq. (|25|) for 
the structure function F2 in charged current vFe scattering at 
a) = 5 GeV^ and b) = 20 GeV^ The sohd curve shows 
the result of our analysis of NuTeV data (one representative 
of 'fit A2') divided by the results obtained with the Base-1 
free-proton PDFs; the uncertainty from the A2 fit is repre- 
sented by the yellow band. For comparison we show the cor- 
rection factor from the Kulagin-Petti model [ij (dashed-dot 
line), HKN07 Q (dashed-dotted line), and the SLAC/NMC 
parametrization (dashed line) [l^] . 



D. Correction Factors for F^{x,Q^) and Fi{x,Q'^) 

We now compute the nuclear correction factors for 
charged current neutrino-iron scattering. The results for 
v-Fe are shown in Fig. [51 and those of u-Fe are shown 
in Fig. [9] The numerator in Eq. (|25|) has been computed 
using the nuclear PDF from fit 'A2', and for the denomi- 
nator we have used the Base-1 PDFs. For comparison 
we also show the correction factor from the Kulagin- 
Petti model [13 (dashed-dotted), and the SLAC/NMC 
curve (dashed) [IjI which has been obtained from an 
A and Q^-independent parameterization of calcium and 
iron charged-lepton DIS data. 



Due to the neutron excess in iron,^^ both our curves 
and the KP curves differ when comparing scattering 
for neutrinos (Fig. [8|) and anti- neutrinos (Fig. [9|); the 
SLAC/NMC parameterization is the same in both fig- 
ures. For our results (solid lines), the difference between 



Note that the correction factors shown in Figs . [Sl and l9l arc valid 
for the case in which the data have not been corrected for the 
neutron excess in iron. For data that aheady have been corrected 
for the neutron excess one should, for consistency, compute the 
/?- factors using A = 56, Z = 28 in equation Eq. l[2}. The mag- 
nitude of the difference between the i?-factors in these two cases 
{Z = 26 vs. Z = 28) is typically a few percent. 
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FIG. 9: The same as in Fig. |8] for vFe scattering. 



the neutrino and anti-neutrino results is relatively small, 
of order 3% at a; = 0.6. Conversely, for the KP model 
(dashed-dotted hues) the u~D difference reaches 10% at 
X ~ 0.7, and remains sizable at lower values of x. 

To demonstrate the dependence of the R factor on the 
kinematic variables, in Figs. [8] and Fig. [9] we have plot- 
ted the nuclear correction factor for two separate values 
of . Again, our curves and the KP model yield dif- 
ferent results for different values, in contrast to the 
SLAC/NMC parameterization. 

Comparing the nuclear correction factors for the F2 
structure function (Figs. [8] and Fig. [9|) with those ob- 
tained for the differential cross section (Fig. [7|, we see 
these arc quite different, particularly at small x. Again, 
this is because the cross section (Pa is comprised of a dif- 
ferent combination of PDFs than the F2 structure func- 
tion. In general, our i?-values for F2 lie below those of 
the corresponding i?-values for the cross section da at 
small X. Since da is a linear combination of F2 and i^3, 
the i?-values for F3 (not shown) therefore lie above those 
of F2 and da. Again, we emphasize that it is important 
to use an appropriate nuclear correction factor which is 
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FIG. 10: Predictions (solid and dashed line) for the structure 
function ratio F^^/F^ using the iron PDFs extracted from 
fits to NuTeV neutrino and anti-neutrino data (fit 'A2'). The 
SLAC/NMC parameterization is shown with the dot-dashed 
line. The structure function F^ in the denominator has been 
computed using either the Base-2 (solid line) or the Base-1 
(dashed line) PDFs. A nuclear correction factor for deuterium 
has been included in the Base-1 calculation [lol |. 



matched to the particular observable. 

As we observed in the previous section, our results have 
general features in common with the KP model and the 
SLAC/NMC parameterization, but the magnitude of the 
effects and the x-region where they apply are quite dif- 
ferent. Our results are noticeably flatter than the KP 
and SLAC/NMC curves, especially at moderate- a; where 
the differences are significant. The general trend we see 
when examining these nuclear correction factors is that 
the anti-shadowing region is shifted to smaller x values 
and any turn-over at low x is minimal given the PDF 
uncertainties. In general, these plots suggest that the 
size of the nuclear corrections extracted from the NuTeV 
data are smaller than those obtained from charged lepton 
scattering (SLAC/NMC) or from the set of data used in 
the KP model. We will investigate this difference further 
in the following section. 



E. Predictions for Charged-Lepton Fi''^'/F^ from 
iron PDFs 

Since the SLAC/NMC parameterization was fit to 
P-i'^ /PF foi' charged- lepton DIS data, we can perform 
a more balanced comparison by using our iron PDFs to 
compute this same quantity. The results arc shown in 
Fig. [10] where we have used our iron PDFs to compute 
F/^"^, and the Base-1 and Base-2 PDFs to compute F^ . 

As with the nuclear correction factor results of the pre- 
vious section, we find our results have some gross features 
in common while on a more refined level the magnitude of 
the nuclear corrections extracted from the CC iron data 
differs from the charged lepton data. In particular, we 
note that the so-called "anti-shadowing" enhancement at 
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X [0.06 — 0.3] is not reproduced by the charged current 
(anti-) neutrino data. Examining our results among all 
the various R[0] calculations, we generally find that any 
nuclear enhancement in the small a; region is reduced 
and shifted to a lower x range as compared with the 
SLAC/NMC parameterization. In fact, this behavior is 
expected given the comparisons of Figs. [TH3] which show 
that at a; ~ 0.1 the cross sections obtained with the base 
PDFs are not smaller than the 'A' and 'A2' fitted cross 
sections. Furthermore, in the limit of large x (x ^ 0.6) 
our results are slightly higher than the data, including the 
very precise SLAC-E139 points; however, the large theo- 
retical uncertainties on in this x-region (see Fig. [5]) 
make it difficult to extract firm conclusions. 

This discussion raises the more general question as to 
whether the charged current (y-Fe) and neutral cur- 
rei it (1 ^-Fe) correction factors are entirely compatible 
d, 0, m, m, m, [5^ . There is a priori no requirement 
that these be equal; in fact, given that the v-Fe pro- 
cess involves the exchange of a and the I'^-Fe process 
involves the exchange of a 7 we necessarily expect this 
will lead to differences at some level. To say definitively 
how much of this difference is due to this effect and how 
much is due to the uncertainty of our nuclear PDFs re- 
quires further study; in particular, it would be interest- 
ing to extend the global analysis of nuclear PDFs to in- 
clude neutral current charged-lepton as well as additional 
charged current neutrino data. Here, the analysis of ad- 
ditional data sets such as the ones from the CHORUS 
experiment [13, H^l (neutrino-lead interactions) should 
help clarify these questions. We are in the processes of 
adding additional nuclear data sets to our analysis; how- 
ever, this increased precision comes at the expense of 
introducing the "A" degree of freedom into the fit. 

V. CONCLUSIONS 

We have presented a detailed analysis of the high 
statistics NuTeV neutrino-iron data in the framework 
of the parton model at next-to-leading order QCD. 
This investigation takes a new approach to this prob- 



lem by studying a single nuclear target (iron) so that we 
avoid the difficulty of having to assume a nuclear "A"- 
dependence. In this context, we have extracted a set 
of iron PDFs which are free of any nuclear model depen- 
dence. By comparing these iron PDFs with "free proton" 
PDFs, we can construct the associated nuclear correction 
factor R for any chosen observable in any given {.t, Q^} 
kinematic range. 

While the nuclear corrections extracted from charged 
current v-Fe scattering have similar characteristics as 
the neutral current l^-Fe charged-lepton results, the de- 
tailed X and behavior is quite different. These results 
raise the deeper question as to whether the charged cur- 
rent and neutral current correction factors may be sub- 
stantially different. A combined analysis of neutrino and 
charged-lepton data sets, for which the present study pro- 
vides a foundation, will shed more light on these issues. 
Resolving these questions is essential if wc are to reliably 
use the plethora of nuclear data to obtaining free-proton 
PDFs. 
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